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Abstract. Analytical estimates of the mass and radial dependence of 
the rates of galaxy mergers and of tidal interactions are derived for clus- 
ters and groups of galaxies, taking into account the tides from the system 
potential that limit the sizes of galaxies. Only high mass galaxies undergo 
significant major merging before being themselves cannibalized by more 
massive galaxies. Strong tides from the group/cluster potential severely 
limit the merger /tide cross-sections in the central regions, and while tides 
are most efficient at the periphery, one should see merging encounters fur- 
ther inside rich clusters. 



1. Introduction 

Mergers of galaxies in slow collisions and tidal interactions in rapid collisions 
are two key dynamical processes that occur in groups and clusters of galaxies. 
Cosmological A^-body simulations are beginning to approach the resolution nec- 
essary to study galaxy dynamics in groups and clusters (see Moore, in these 
proceedings). Moreover, mergers and tidal collisions leave significant observa- 
tional signatures, in the form of tidal tails, asymmetries and generally disturbed 
morphologies and internal kinematics (see Amram, in these proceedings). Also, 
galaxy merging is an essential mechanism for driving elliptical galaxy morpholo- 
gies given disk-like progenitors. As such, an understanding of galaxy merging is 
very important for semi-analytical modeling of galaxy formation. 

In this review, I compute analytically the rates at which a galaxy of given 
mass and position within a cluster or group with a Navarro, Frenk, fc Whit^ 



(|T99|, NFW) potential under goes slow major mergers with lower iiass galax- 



ies and rapid tidal encounters. Since the collision cross-sections are strongly 
modulated by the tides from the group/cluster potential, I begin with a simple 
formalism for estimating the tides from the system potential. Note that I will 
not consider ram pressure stripping on galaxies. 



2. Tides from the cluster/group potential 



The potential of a cluster can exert a strong differential force on a galaxy orbiting 
within it, but these tides are strongly dependent on the galaxy orbit. 
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A galaxy on a nearly circular orbit is likely to be tidally locked, as the Moon 
is with respect to the Earth. In this case, the tidal force is simply ( |King, 1962 ) 



tide 



A 



GM{R) 



(1) 



and equating Ftide in eq. (0) to the force, / = Gm{r)/r'^ that a galaxy exerts on 
one of its stars yields for r <^ R a galaxy tidal radius rt that satisfies a velocity 
modulated density criterion (see Mamon, 1995| ): 



Pgin) = Pcl(-R) 



+ 



(2) 



where Vp and Vr 



are respectively the galaxy's velocity at pericenter and the 
cluster's circular velocity. The term in brackets in eq. (§) is 2 for circular orbits. 
Moreover, for a singular isothermal law p ~ for both galaxy and cluster, one 
finds (see White, 1983| ) that for circular orbits the galaxy size is proportional 
to its clustocentric radius: rt = 2~-^/^(fcirc/K;irc)-R) where Vdrc is the circular 



velocity of the galaxy. For general density profiles, writing p{r) 



one obtains 



n/R 
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Merritt (1984 ) used a similar circular-tide criterion to argue that galaxies are 
strongly tidally limited by the cluster potential. 

However, cosmological infall imposes elongated orbits. A galaxy on an 
elongated orbit experiences a strong tide during its rapid, hence short, passage 
at pericenter, prompting Ostriker, Spitzer, fc Chevalier (1972| ) to introduce the 
term tidal shock. During this shock, a star in the galaxy experiences a velocity 
impulse 



GM{Rp) r R 



i?3 



GM{Rp) 



(4) 



where we neglected the centrifugal term in Ftide, becaus e the galaxy f alls in too 
fast to be phase locked. A more precise calculation by Spitzer (1958 ), who in- 
troduced the impulsive approximation where the point-mass perturber moves at 
constant V, produces the same relation as in eq. (Q) with a constant of order 
unity. The impulse approximation can also be applied to extended perturbers 
(lAguilar fc White, 1995| ; ^ amon, 1987| ). Recently, [ Gnedin, Hernquist, fc Osl 
triker (1999]) applied the impulsive approximation to the more realistic Hernquist 
(19901 ) potential and found a dependence of Av matching that of eq. with 
a very small dependence on Rp and again the constant is found to be of order 
unity for elongated orbits (with order of unity changes when they performed 
orbit-integrated — instead of straight-line — tidal calculations). 

The tidal radius can then be defined as that where the energy increment 
caused by the tidal perturbation is equal to the binding energy ( [White, 19831) . 
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With E ~ Gm{rt)/rt and AS ~ (At>)V2, one obtains another velocity modu- 
lated density criterion 



Pgin) ^ Pc\{Rp) 

which for any density profile yields 
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Figure 1. Velocity modulation of normalized tidal radii for given 
orbit pericenter. The solid curve shows impulsive tides for elongated 
orbits (eg. p). The dash ed curves show circular tides (eq. j^) in a 
Navarro, Frenk, &: White potential for 0.1 [upper curve) and 1 [lower 
curve) scale radii. The dotted curve shows the circular tide within a 
homogeneous core. 

Figure |l| shows the effects of velocity modulation on the tidal radii. The 
effective tidal radius should be taken as the largest of the circular and impulsive 
regimes (otherwise one would be left with discontinuities in the transition Vp 
from near-circular to elongated orbits). Hence, with homogeneous cores, circu- 
lar tidal theory produces increasingly smaller tidal radii for orbits of increasing 
but low elongation, as Merritt &: White (1987) found in their A^-body simula- 



tions. With cuspy cores as in the NFW profile, low-elongation non-circular orbits 
experience tidal shocks instead. The orbit elongation, Rp/Ra — 0.2, found in 
|Ghigna et al. 's (1998) high-resolution cosmological simulations corresponds to 
Vp/Vcirc[Rp) ^ 1-5 - 2.7, roughly yielding 



n 



(1.2 - 1.5) Rp 



[1.2 - 1.5) Rp{^ 

\Vcl 



(7) 
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as for singular isothermal models, where Vg and Vd are the mean galaxy and 
cluster velocity dispersions, respectively. If ~ Vcirc{Rp) and if galaxy and 
cluster density profiles are self-similar, then tides from the cluster potential 
would force the simple relation m{rt)/m{r^i^) = M {Rp) / M {R^\^) . In fact, using 
yp/yc\rc{Rp) expected for the NFW profile, assuming (i?) ~ 4fi n^and adopting 
the departures from self-similarity in the NFW profiles noted by ^ 



avarro, Frenk, 



& White (1997, lower mass NFW profiles are more centrally concentrated), one 
finds 



m{rt) 



M{R) 



M{R, 



R 



(8) 



where for clusters {vd = 1000 km s^ ) and groups {vd = 300kms~ ) we respec- 
tively have bm = 0.50 and 0.82, Ur = 0.58 and 0.52, and br = 0.57 and 0.78 
(eq. 1^ is accurate to better than 5% for R > 0.05i?vir)- 



3. Galaxy merger rates in clusters and groups 



3.1. Global merger rates for equal mass galaxies 

The rate of mergers is obtained by integrating over velocities the merger cross- 
sections: 

1 (PN f°° 
^ ^ ^ dtdV " ^^'^(^^^ " Jq dvf{v)vs{v) , (9) 

where s{v) = 7r[pcrit(f )]'^ is the merger cross-section and f{v) is the distribution 
of relative velocities (with f{v)dv = 1). Hence, nk = dN/dt is the rate at 
which a galaxy suffers a merger. Within the virialized regions of clusters with 
ID velocity dispersion Vd, the velocity distribution is a gaussian with standard 



deviation 2-^/'^Vd: f{v) = 2 -'^tt-'^/'^v~^^v^ exp[-v ^ / (Av-^^)]. 



Roos fc Norman (19V 9 '), [Aarseth &: Fall (1980D and [Farouki fc Shapiro (T982|) 
have established merger cross-sections from very small A^-body simulations of 
galaxy collisions, that were based upon the parameters at closest approach. 
The maximum distance of closest approach, r™*^^, was 4 ( [Aarseth fcTallD or 11 
(Farouki & Shapiro) times the mean galaxy half- mass radius, r^. Mamon (1992 ) 



used the Roos &: Norman cross-section with the Aarseth Sz Fall scaling to derive 
a merger rate. 

However, the gaussian approximation for the relative velocity distribution 
implies that the cross-sections used in eq. (| ) are based upon impact parameters 



(at infinity) and not at closest approach. [Makino &: Hut (1997| ) derived merger 
cross-sections using high resolution A^-body simulations of colliding galaxies with 
more realistic density profiles. Their cross-sections are expressed in terms of 
impact parameters. 

Krivitsky & Kontorovich (1997) used a simple gravitational focusing recipe 
to connect the cross-section at closest approach (which they assumed to be 
independent of pericentric velocity), in units of the velocity at infinity, to the 
cross-section at infinity (without making any a ssumption on th e potential energy 
of interaction of the colliding pair) . Note that Makino &: Hut| show that r™^^ > 
lOr/j, while pTrivitsky &: Kontorovich] argue that it is approximately the sum of 
the galaxy radii (which are ill-defined). 
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where is the 



Figure ^ shows the dimensionless merger rates k/{rf^Vci, 

mean galaxy internal velocity dispersion, derived from |Mamon (199^ ), |Makino & 
Hut (1997 ), and Krivitsky Sz Kontorovich (19971 ) • I rescale the rates of [Krivitsky 
fc KontorovichI in terms of half-mass radii using R = 9rh (to obtain merger rates 



similar to those of Makino &: Hut] ), i.e., r'^^^/rh = 18 
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Figure 2. Dimensionless merger rates k/{r\vg) (see eq. (P) as a func- 
tion of the ratio of cluster to galaxy velocity dispersion. The dotted, 



solid and dashed curves respectively correspond to the rates of Mamon 
(19921), [Makino fc Hu^ (|1997| using [Hernquist, 1990| model galaxies) 



and Krivitsky & Kontorovich (1997, rescaled vertically, because they 
reason in terms of galaxy radii rather than galaxy half- mass radii). 



The agreement between the merger rates of Krivitsky &: Kontorovich and 
Makino fc Hut| is remarkable, given the very simple analytical formulation of the 
former authors (but again this required a rescaling, or, in other words, a choice 



of r^^^/vfi = 18). The very low merger rate of Mamon (1992 ) in the group 
regime {vc\ 
model. 



Vg) is a consequence of the lack of gravitational focusing in that 



The merger rates of Mamon (1992) and Makino & Hut agree to within 
15% in the cluster regime (fci ^ 4ug). This agreement is almost fortuitous since 
Mamon shows that the merger rate with the Roos & Norman cross-section scales 
as (r™^^/r/i)^, while this ratio is very different in Makino & Hut's cross-section. 
In any event, in the cluster regime the merger rate can then be written 



v: 



cl 



(10) 



where a ~ 8 ( Mamon, 1992| ). With Sv^ O.AGm /rh ( ^pitzer, 1969| ), appro- 
priate for the Hernquist model, the Makino fc Hut rate translates to a = 12. 
Figure |2| shows that k ~ w^^, whichever merger cross-section is used. In fact, 
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it is easy to show that for any merger cross-section rapidly decreasing with in- 
creasing velocity, the merger rate should scale as v~,^ for Vc\ ^ Vq, as first found 



by Mamon (1992 ) for the Roos &: Norman] cross-section. 



The important conclusion of Figure is that for given galaxy parameters, 
the merger rate is roughly 100 times lower in rich clusters than in poor groups 
of galaxies. 

3.2. Merger rates for difTerent masses 

If the critical merging velocity fcrit is a function of rp/r^ ( [Aarseth &: Fall, 198C ; 



Farouki &: Shapiro, 1982| ), it is easy to show that k ^ r\, and if t^crit is a function 



of rp/{rh), then k ~ (see Mamon, 1992). Similarly, it is reasonable to 

expect that k ~ {Vg)'^. Then, given eg. (|1[1D and that m ~ r|, the rate of 
mergers of a galaxy of mass m with a galaxy of mass Am 

A given galaxy undergoes mergers with other galaxies at a rate 

_ /"Amax 

TZ = nk{m) = / k{m, Am) n{\m) d{\m) , (12) 

where for major mergers with smaller galaxies (that transform disk galaxies into 
ellipticals), Amin — 1/3 and Amax = !> while for des truction by mer gers with 
larger galaxies, Amin = 1 and Amax — > oo. Adopting a ^chechter (1976| ) form for 
the mass function of galaxies, n[m) = {n^,/m^) x~" exp(— x), where x = m/m^,, 
egs. (InD and (|l|) yield 



aG'^n^ml 
16^3 



n = nk= 3 * K{m/m,) , (13) 



6 

i^major(a;) = x3-°^Min(j,7-j) [r(l + j73- a, x/3) - r(l +i/3 - o, x)] , (14) 

j=0 

6 

i^destr(x) = x=^-°^Min(j,7-i)r(l+i/3-a,x) . (15) 

j=0 

Figure ^ shows the expected number of major mergers with smaller galaxies 
and destruction by mergers with larger galaxies that a galaxy of a given mass 
should expect in a Hubble time if it sits in a typical location of a rich cluster, 
assuming a constant rate in time. The rise in merger rates at low mass reflects 
the rise of merger cross-section with mass, while the decrease at high mass is 
caused by the sharp decrease in the galaxy mass function yielding few galaxies to 
merge with. Figure ^ clearly indicates that the probability of merger for a given 
galaxy is always small. Moreover, low and intermediate mass galaxies (m < m*j 
are usually cannibalized before undergoing major mergers. 
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Figure 3. Number of major mergers with lower mass galaxies 
(eqs. |13| and [Q, thick curves) and higher mass galaxies (eqs. and 
||l5|| , thin curves) extrapolated to one Hubble time versus galaxy mass. 
Dotted, solid and dashed curves are for a = —1.1, —1.3 and —1.5, 
respectively. The normalization assumes a = 12, Vd = 1000 km s~^, 
n, = 200nf with n^°'^ = 0.013 /i^Mpc'^ QMarzke et al, 19981 ), and 
m* = 0.1 (M/L)ci 4 = 3 X lO"/i-iM0. 



3.3. Variation of merger rates with position in cluster 

One can go one step further and predict the variation of the merger rate with 
position in the cluster. From eq. (p!3|), the merger rate scales with radius as 



aG^m^^^{R)pciiR) 

T^iR,m) = T— prr -^—K{m m^) , (16) 

16T{2-a,Xm)v^i{R) 

where p,{R) = m{rt) /m{r^-„) (see eq. §), m^, is the mass at the break of the 
field galaxy mass function and Xm is the minimum galaxy mass in units of m*. 

Assuming a mass density profile p ~ R^^ and arguing that cluster galaxies 
are severely tidally truncated by the cluster potential as rgai ~ R (i.e. with 
eq. and assuming constant Rp/Ra), [Mamon (1992 ) showed that if galaxies 



also follow the law p ~ r~^, then their masses obey m ~ R^~^. Note that this 
sharp scaling of galaxy size with clustocentric distance is now confirmed in high 
resolution cosmological simulations of clusters ( Ghigna et al., 199^ ). Hence, the 



radial variation of merger rates are strongly modulated by potential tides. 

By writing n{R) ~ p{R)/m{R) ~ R^^ , I derived nk ~ hence a 

higher merger rate inside the cluster, with a slope agreeing perfectly with the 
observed elliptical fraction ( [Whitmore, Gilmore, &: Jones, 1993]) , given / ? = 9/4 
as predicted in early models of cluster formation ( |Bertschinger, 1985| ). The 
derivation above has one flaw: although galaxy masses were correctly scaled 
to increase with R, I forgot to scale the fraction of cluster mass lying within 
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galaxies in the same way. Therefore, one really expects n{R) ~ p{R) ^ R ^ and 
nk ^ _r3-3/3/2 yigifiing a slope (iln(nA;)/(ilni? = —3/8 for j3 = 9/4 and a null 
slope for (3 = 2, both in disagreement with the logarithmic gradient of elliptical 



fraction found by Whitmore et al 



One can use the more realistic NFW density profiles to estimate the ra- 
dial dependence of the merger rates. An essential parameter is (Rp/R), which 
measures the effectiveness of the tides from the cluster potential. Because the 



dynamical friction time scales as M/m times the orbital time (Mamon, 1995), 
orhit circularization, which to first order operates on a dynamical friction time 
scale, should be very slow for galaxies falling onto clusters, but fairly effective 
for galaxies falling into small groups. 

Figure ^ shows the predicted number of major mergers in rich clusters and 
small groups, extrapolated over a Hubble time, using the non-self-similarity of 
the NFW profiles, an exact scaling of the typical galaxy mass with radius 
(see eq. ||8|), and partial orbit circularization in groups. 
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Figure 4. Number of major mergers a given galaxy undergoes with 
lower mass galaxies (eq. [jl6|, with eqs. [13| and [Q), extrapolated to 
one Hubble time, versus clustocentric radius in (left) an NFW cluster 
with Vci = 1000 km s~^ and (right) a group with Vd = 300kms~^, where 
the galaxy mass function has a = 1.3. The solid thick and dashed thick 
curves represent the expected number of major mergers for galaxies of 
mass m = m^"^'^ = nopc/[nf^'^T{2 - a, Xm)] = bxlO^'^h'^MQ and m = 
0.1 m^'^^'^ = 5 X W^^h^^MQ, respectively. The solid thin and dashed 
thin curves refer to galaxy masses m = m^{R) and m = 0.1 m^{R), 
respectively. 
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The merger rates for constant mass galaxies fall off sharply at small clusto- 
centric radii, simply because tidal truncation of galaxies is so severe that there 
are no galaxies left that are massive enough to produce a major merger with our 
test galaxy. In general, the merger rates are maximum for intermediate radii 
for given test galaxy masses, and at low radii for fixed m/m^{R) (recall though 
that low mass galaxies get cannibalized before they can undergo a major merger 
with a smaller galaxy). In any event, Figure ^ confirms that mergers are ineffec- 
tive in clusters, hut very effective in small groups. Note that without resorting 
to partial orbit circularization within groups, the expected number of mergers 
in groups is somewhat less than expected from the scaling, because, rela- 
tive to clusters, the more concentrated NFW profiles in groups lead to stronger 
modulation of the merger rate by the potential tides. 



4. Collisional tidal stripping in clusters and groups 

Because non-merging galaxy collisions are by essence rapid, they can be treated 
as tidal shocks, and it is reasonable to assume that for tidal features to be visible, 
one requires Ai? > ^\E\, where 7^1- Hence, > (87) ^^^f^. Denoting p and 
V the separation and relative velocity at pericenter, and fcirc,p(p) the circular 
velocity of the perturbing galaxy out to p, eq. (Q) leads to a critical impact 
parameter 

r v^- 

' ^circ,p (-\^\ 

- (37)1/2 v,V ' ^ 

where we note that Vcirc,p5 is almost independent oip for realistic density profiles 
for the perturbing galaxy. Then, integrating the cross-sections derived from 
eq. ( p!7D , the rate of tidal interactions is 

1/2 / ^ \ %4 

37 \V„ f cl 



and is virtually independent of the test galaxy parameters. Integrating eq. (|l^ 
over perturber mass and remarking that both the galaxy and the perturbers are 
tidally limited by the cluster potential, one obtains using eq. (^) 



- r(7/3-a,x^) (v^,,^\ p^^{R)ij,y^{R) 



n 2 



^irc i^p) 



47rV2-y Gr(2-a,Xm) \ Vg J ni^ pci{Rp) VdiR) 

(19) 

where again fi{R) = m{rt) /m{r^iy-) (see eq. [P), Vdrc,* is the circular velocity 
at the virial radius for a field m^, galaxy. Again, the rate of tidal encounters is 
independent of the galaxy mass. Note that the fi'^^^{R) dependence of the rate of 
tidal encounters illustrates the strong modulation of this rate by the tides from 
the cluster potential. 

Figure ^ shows the expected (eq. [|19|) number of strong tidal collisions for 
galaxies in clusters and groups, with {R/Rp) as in Figure ^. Although groups 
are preferential sites for strong tidal encounters, galaxies in the outskirts of 
clusters should also witness such interactions. However, the signature of tidal 
interactions lasts of order lh~^ Gyr, so that the fraction of galaxies in clusters 
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Figure 5. Number of strong (7 = 1/3) tidal encounters a given galaxy 
undergoes (eq. fl^), extrapolated to one Hubble time, versus clusto- 
centric radius in an NFW cluster with Vc\ = 1000 km s~^ {solid curve) 
and a group with Vd = 300kms~^ {dashed curve). The galaxy mass 
function has a = 1.3 with field = 5 x 10^^ h^^ Mq. 



and groups that are currently undergoing tidal interactions is roughly one-tenth 
of what is displayed in Figure ^. 



5. Discussion 



The strong radial dependence of galaxy masses, predicted by the tidal theory 
(eq. ^), is clear in the cosmological simulations of phigna et al. (1998 ). Should 
we then witness inverse luminosity segregation in clusters where, outside of the 
core, galaxies become more luminous towards the cluster periphery? Indeed, 
Adami, Biviano, fc Mazure (1998| ) found a weak trend of mean galaxy magnitude 



versus radius for an ensemble of clusters, although they worry that this trend is 
caused by observational bias. It may be that incompleteness of the observational 
samples is washing out the trend rather than creating it. 

The lack of mergers in present-day rich clusters has been noted in cosmo- 
logical simulations of clusters (Ghigna et al., 1998). Figure^ shows that in rich 
clusters, mergers are at best marginally probable for high mass galaxies lying in 
the cluster body. Given that high mass galaxies are rare, such merging will be 
difficult (but not impossible) to detect observationally or in simulations. 

From their Ha prism surveys of galaxies in clusters, Moss and co-workers 
( IMoss fc Whittle, 1993| ; |Moss, Whittle, fc Pesce,l99l ; [Bennett fc Moss, 1998|) 
note ~ 30% of spiral galaxies in rich clusters exhibit a compact Hq morphology 
and roughly half of these tend to be morphologically disturbed and have nearby 
neighbors. Another half of these compact Ha emission galaxies are in the cluster 
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core, and presumably those have no tidal companions, suggesting that they are 
harassed by the cluster potential. But the first half, outside the core are probably 
bona fide cases of strong tidal interactions within clusters, leading to an overall 
frequency of 20% of all galaxies outside the cores of clusters. It remains to be 
seen if they are associated to substructures such as infalling groups. Correcting 
Figure ^ for the 1 Gyr duration of these tidal features produces an absolute 
frequency of tidally interacting galaxies in clusters in rough agreement with 
that found by Moss and co-workers. 
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